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Abstract Self-energy corrections involving logarithms of the 
parameter Za can often be derived within a simplified ap- 
proach, avoiding calculational difficulties typical of the prob- 
Ph. lematic non-logarithmic corrections (as customary in bound- 

ed ! state quantum electrodynamics, we denote by Z the nuclear 

charge number, and by a the fine-structure constant). For some 
logarithmic corrections, it is sufficient to consider internal prop- 
erties of the electron characterized by form factors. We provide 
a detailed derivation of related self-energy "potentials" that 
give rise to the logarithmic corrections; these potentials are lo- 
cal in coordinate space. We focus on the double-logarithmic 
two-loop coefficient Bq2 for P states and states with higher an- 
gular momenta in hydrogenlike systems. We complement the 
discussion by a systematic derivation of Bq2 based on nonrela- 
tivistic quantum electrodynamics (NRQED). In particular, we 
find that an additional double logarithm generated by the loop- 
after-loop diagram cancels when the entire gauge-invariant set 
of two-loop self-energy diagrams is considered. This double log- 
arithm is not contained in the effective-potential approach. 
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1 Introduction 



Lamb-shift measurements and related theoretical calculations for bound atomic systems with in- 
creasing accuracy have historically provided accurate tests of quantum electrodynamics (QED), 
and the measurements have recently been improved in accuracy beyond previous limits 
In order to account for a theoretical description, corrections of various physical origin (one-loop 
self-energy and vacuum polarization, two-loop, and higher order radiative, recoil, radiative-recoil, 
nuclear-size corrections) have to be evaluated Q. 

Here, we focus on logarithmic self-energy corrections which are evaluated within the Za- 
expansion Q. Within the analytic treatment, self-energy radiative corrections can be taken into 
account by means of a nonanalytic expansions in powers of the fine-structure constant a, the 
product of Za and the logarithm ln[(Za)~ 2 ] (Z is the nuclear charge number). The expansion 
in powers of a corresponds to the loop-expansion in the framework of the usual perturbative 
treatment for QED. The higher order terms in powers of Za and ln[(Za) -2 ] are related to 
atomic-physics effects; they are referred to as the "binding corrections". 

The purpose of this investigation is twofold: first, to illustrate how Lamb-shift "potentials" that 
give rise to the logarithmic corrections can be derived within the context of bound-state QED, 
and second, to provide a rigorous and detailed derivation of the B§2 double-logarithmic two- 
loop self-energy coefficient for P states and states with higher angular momenta. The P-state 
coefficient Bq2 has already appeared in the literature || ; however the derivation has been rather 
sketchy. 



2 Modified Dirac Hamiltonian, One— Loop Corrections and A41 

It has been observed by many authors (e.g. 0^9|) that a rather important class of self-energy 
radiative effects for bound states can be described by a modified Dirac Hamiltonian (K = c = 
eo = 1), 

= a[p-eFi(A)A]+^m + eFi(A)</> 

+F 2 (A)^(ij-E-l3a-B), (1) 

which approximately describes an electron subject to an external scalar potential <f> = cp(r) and 
an external vector potential A = A(r) (the vector potential vanishes for a point nucleus that 
gives rise to a static Coulomb potential; we may neglect the nuclear magnetic field and the 
hyperfine structure). We have 

e <p(r) = eA Q (r) = -— (2) 
r 

in coordinate space, which corresponds to 4>{q 2 ) = —AirZa/q 2 in momentum space. In this 
article, following the commonly accepted convention, the function (p(r) and its Fourier transform 
<f){q 2 ) are denoted by the same symbol 4>. We avoid possible ambiguities by denoting with r and 
r the arguments in coordinate space and with q or p those in momentum space. The argument 
A = (d/dr) 2 of the electron form factor F\ in Eq. ([!]) is to be interpreted as a Laplacian operator 
acting on all quantities to the right, but not on the wave function of the bound electron ip(r). 

Equation (jlj) entails a replacement of the binding Coulomb potential as 

e<f){r) -> eFi(A)<£(r) 
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and leads to a correction to the Coulomb potential AVc(r) according to 



AI r (D- i.f,.(A)- I] ( 



in coordinate space, and 



AF C (<7 2 ) = [iM-Q 2 )) - 1] 



4-7T 



(3) 



(4) 



in momentum space. In first-order perturbation theory, this gives rise to the following pertur- 
bative correction which we write down in coordinate and momentum space, 



AEx 



\AV c (r) | V) = (i/> | [-Pi (A) - 1] e 1 V) 



d d r V 7 (r) 



[Fi(A)-l] 



d 3 p /" d 3 p' 



(2vr) 3 7 (2tt) 



^+(p') 



-4vrZa\ 



(5) 



with q = p' — p. An expansion of the electron form factor Fi in terms of its argument gives rise 
to higher-order terms in the Za-expansion, because the atomic momentum is of the order of Za 
in natural units. Therefore, within the Za-expansion, it is admissible to expand both the bound- 
state Dirac wavefunctions tp in powers of Za (the leading-order term is then the Schrodinger 
wavefunction) , as well as the electron form factor in Eq. (||) in powers of its argument. 

The one-loop (1L) self-energy correction for S states within the Za-expansion reads 



AE, 



(1L) 



SE 



-) (Za) 4 ™ (A 4 i H(Za)- 2 ] + A m + K) 



(6) 



,7r/ n° 

where the remainder 1Z is of the order of O(Za), m is the electron mass and n is the principal 
quantum number. 

As indicated in Eq. (||), the form factor -Fi(A) in momentum space assumes arguments according 



to the replacement A 



-(p'—p) in momentum space. With the convention q = q^q^ 



(q ) —q , the evaluation of the radiative corrections to the binding Coulomb field is mediated by 
space-like virtual photons (q° = 0), and the momentum transfer can be written as: q 2 = —q 2 = t 
(this is consistent with the conventions employed in 



W) 



The form factor F\(t) can be expanded in powers of a, which corresponds to the loop expansion. 
According to Eqs. (1.2) and (1.20) of [pLCj ] , we have up to two-loop order: 



a\2 

7T 



F^(t) + 0(a 3 ) 



with 



B(t) ln-+f[ 2) (t). 
m 1 



1 



A 



A 



Ef\t) In" — + In — B(t) T {2) (t) + Ff } (t) 



r(4) ( 



(7) 



(8) 



(9) 



B(t) [see Eq. (1.18a) of [10|] reads as follows 



where the T are infrared finite (i.e. finite in the limit A — > 0), and the definition of the function 

(10) 



B(t) 



1 + 



t-2m 2 



[1 - Am 2 /tf/ 2 - 1 
t(l - 4m 2 /0 1/2 " (1 - 4m 2 /«) 1/2 + 1 



In 



3r?r 



+ 0{t z 
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In Eq. (||), A denotes the fictitious photon mass. How should the problem of the infrared diver- 
gence of the form factors be interpreted in the context of bound-state QED? The free electron 
can emit an infinite number of infrared photons, because it may undergo transitions between 
free states with infinitesimal energy differences. However, this is not the case for a bound elec- 
tron which has a discrete bound-state spectrum; energy levels are separated from each other by 
intervals of the order of (Za) 2 m (the energy level differences are determined by Schrodinger 
theory). This leads to an infrared cutoff in bound-state QED of the order of A ~ (Za) 2 m. 
Therefore, we may replace A — » (Za) 2 m for the determination of leading logarithms of the 
Lamb shift. At some risk to over-simplification, one may therefore argue that the infrared catas- 
trophe is avoided in a natural way for bound states. For the description of bound states, we have 
m(A/m) ~ — ln[(Za)~ 2 ] within logarithmic accuracy, i.e. neglecting non-logarithmic contribu- 
tions which are given e.g. by A40 coefficients [see Eq. (||)]. 

The focus of the current article is on double-logarithmic corrections which are present from the 
first term on the right-hand side of Eq. (S) . Note that single- logarithmic two-loop corrections are 
not being considered in this article. Corrections of this latter type are generated, for example, 
by the second term on the right-hand side of Eq. (|9|) . 

At this point, it may be helpful to point out that the cutoff of the infrared divergence of QED 
at the "bound-state photon mass" A — > (Za) 2 m is consistent with the matching procedure 
that involves an explicit infrared cutoff e which can be interpreted as an infrared cutoff for the 
bremsstrahlung spectrum |7],§,12]. The procedure is described in some detail in Eqs. (32) - (34) 



of This matching procedure offers an alternative interpretation for the infrared catastrophe: 
the infrared divergence crucially relies on transitions between asymptotically free electron states. 
Any infinitesimally small additional interaction of the electrons within that interferes with the 
emission of bremsstrahlung will avoid the infrared catastrophe and provide an infrared cutoff 
whose order-of-magnitude is determined by the energy scale of the additional external field. 

In combining the result @ with the expansion of B(t) in powers of t, we reproduce the well- 
known expression 



F?\t) 



t 



3m 2 



, A 1 
In h - 

m 8 



+ 0(t 2 



(11) 



Together with the definition of the modified Coulomb potential in Eq. (Q) and the bound-state 
"infrared-cutoff prescription" A — > (Za) 2 m, this leads to the following one-loop (1L) self-energy 
potential 



a 

IT 



3777/ 



M(Za)- 2 ]) 



4irZa\ 



4a 
3m 2 



(Za) ln[(Za) 



(12) 



in momentum space; this translates into a potential 



Av£ L \r) 



— (Za) In (Za) ] — 

6TT m z 



(13) 



in coordinate space. This potential can also be found as Eq. (2) of 13], given there without 
derivation. The first-order one-loop perturbation, evaluated according to Eq. (||), reads 



AE[ 1L) = (^|A^ 1L) (r)|V) 



4a 
3^ 



(Za) A ^ln[(Za)- 2 ]5 l0 . 



(14) 



This correction is nonvanishing only for S states (I = 0) , and it reproduces the leading logarithmic 
A41 coefficient as given in Eq. (||). It may be interesting to point out that since \ip(r = 0)\ 2 = 
(Za) 3 (m 3 /tt) 5iq, where m r is the reduced mass of the system, the correction (14) also has the 



4 



correct reduced-mass dependence (this is of relevance for systems like positronium and pionium). 
In the limit of a large nuclear mass, we have of course m = m r . 

Note that the potential (13) is local in coordinate space. In contrast, the nonrelativistic (NR) 
one-loop self-energy operator (as well as its relativistic counterpart which assumes a slightly 
more complicated form) may be expressed in the length-gauge form as [cf. Eq. (29) of [[T|])], 



,(1L) 
J NR 



r,r') 



2a 
3^ 



du uj 3 r' 



1 



H-E + co 



(15) 



where e is the upper cutoff for the photon energy originally introduced in p2|] . The self-energy 
operator (|l5|) involves two spatial coordinates. The locality of the potential expresses the 
fact that the high-energy virtual photons which mediate the form-factor corrections in Eq. (|l|) 
act on a relativistic length scale given by the Compton wavelength of the electron which is 
smaller by one order of Za than the atomic length scale given by the Bohr radius. 



3 Effective Local Potential for 
Two— Loop Corrections and Bq2 



In combining the result (|9]) with the expansion of B(t) in powers of t [see Eq. (|1C|)] and the 
modified Coulomb potential in Eq. (B), and using the bound-state "infrared-cutoff prescription" 
A — > (Za) 2 m, we obtain the following two-loop (2L) self-energy potential 



AV^<~ 2 



id 2 ) 



a\2 1 



7T 



Q\2 1 



3 m' 



- — ln 2 [(Za) 
ttJ 18 iy ' 



ln 2 [(Za)~ 2 ] 
o, 47r Za , 



Air Za 



mr 



(16) 



This correction has previously appeared as Eq. (3) of 
Fourier transformation, we have 



without a detailed derivation. After 



Avg h \r) 



9 W 



ln 2 [(Zft)" 



7T AS^Hr) 



mr 



(17) 



which is a highly singular potential in coordinate space. Its expectation value on S states diverges, 
giving rise to a further logarithm, and we will not discuss here the associated problems, which 
have recently attracted remarkable attention [15-22]. 

The first-order perturbation, evaluated according to Eq. (||), reads 

-ft\2 2vrZft 

-7T/ 



AE- 



(2L) 



mr 



In 2 [{Za)- 2 ] A 



)n,l=l,m 



r=0 



(18) 



In Eq. (18), the Laplacian operator acts on a Schrodinger P wavefunction. The following analytic 
result 



A 



11 



1 



m 



(19) 



r=0 3n 

where n is the principal quantum number, has previously appeared in the literature (e.g. |],||]). 
Within the current investigation, we would like to present a complete derivation of the analytic 
expression for this matrix element in the Appendix |A[ Finally, we rewrite the energy correction 
in the form 



AE 



(2L) 



(Za) e m 



nr 



i» 2 [(z*r 2 } yi 



4 n l 



1 



n- 



(20) 
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This double-logarithmic correction originates solely from the two-loop F\ form factor of the 
electron. This corresponds to the diagrams in Fig. |l| (a) and (b). To complete the gauge- invariant 
set, the loop-after-loop diagram in Fig. |] (c) should also be taken into consideration. 

The diagram in Fig. |l| (c) gives rise to a "second-order perturbation" involving to one-loop self 
energies as first-order perturbations (the "irreducible part" of the diagram), supplemented by 
a further term involving the derivative of the bound electron's Green function (the "reducible 
part"). The correction is known to read (see e.g. p3| ) 



4 1L) (£) I £ 



W) w 

E — E,.i 



4 1L) (^) 



1>) , (21) 



where S^ 1 "^ (E) is the renormalized relativistic one-loop self-energy operator, and E is the energy 
of the electron in the state Within the effective-potential approach, the one-loop potential 
( |T3| ) describes the two one-loop self-energy insertions in the first term of (^)- The potential ( |T3| ) 
involves a Dirac delta-function in coordinate space that vanishes on P states, and consequently 
it can be argued that no further double- logarithmic corrections originate from this term ( but see 
the discussion in Sees. |I| and ||). 

The second term in (|2l|), which involves the derivative of the self-energy operator with respect 



to its argument [see also Eq. (2.6) of JT2| or Eq. (2) of ||] and constitutes the reducible part of 
the diagram in Fig. [l| (c), does not give rise to any further double logarithm, either. The first 

factor ftp EjJ L ^(^) tp^j does not create any logarithm for P states in the order of a (Za) 4 '. The 
second factor, which contains the derivative of the self-energy operator, is not separately gauge 
invariant, and consequently, there exists no "effective potential" which could be inserted for 
this term. This is in itself a rather unsatisfactory situation for the effective-potential approach. 
However, it is possible to analyze the logarithm which is generated by the nonrelativistic photon 
integration region in this term. Consider the nonrelativistic "velocity-gauge" form of (|1]) and 
differentiate with respect to the energy, 



±^(e) 



NR 




v_ 

m 



1 



H-E + lo 



V_ 

m 



(22) 



where <f) is the nonrelativistic (Schrodinger) wavefunction, There is only a single logarithm 
ln[e/(Za) 2 m] generated in the integration region lo G [(Za) 2 m, e] which may be extracted by 
replacing 1/(H — E + uj) — > l/u>. The logarithmic term is proportional to the matrix element 
(4>\(p 2 1 m 2 )\(j)) , which is finite on P states. Consequently, no further double logarithms arise from 
the second term of (Ell). 



The two-loop effect for P states is usually characterized by the following semi-analytic expansion 
in powers of Za [cf. Eq. @], 



AE, 



(2L) 
SE 



7T 



in 



- (Zaf [B m + {Za) 2 [B w In 2 (Za) 



-2 



+B 61 \n{Za)- 2 + B m }+K) 



(23) 

where the remainder 1Z is order O(Za) 3 . Using Eq. (|20|) , one can immediately read off the 
two-loop double logarithmic spin-independent coefficient 



Be2(n,l 



i 



27 n 2 

We confirm the result obtained for this correction in Ref. 



(24) 
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00 

Figure 1: The crossed (a), rainbow (b) and the loop-after-loop diagram 
(c) which contribute to the two-loop self-energy for a bound electron. 
The propagator of the bound electron is denoted by a double line. 



4 Double Logarithms and the Loop— After— Loop Diagram 



In the previous section, we have seen that within the effective-potential approach, no double 
logarithm originates in the order (Za) e from the loop-after- loop diagram in Fig. [l] (c). This is 
because, within this approach, we insert the delta-like local potential ( |l~3|) for the two one-loop 
self-energies in the first term of (|2l|). 

However, if we consider the diagram in Fig. |l] (c) within the Coulomb gauge and formulate 
the contribution due to low-energy virtual photons, then we obtain for the irreducible part the 
expression 



A£ L al 



■ > n,l,m 



,(1L) 



1 



H-E 



J NR 



■ > n,l,m 



(25) 



where the nonrelativistic self-energy operator is given by Eq. (jT5|) , and (f) nt i t m is the Schrodinger 
P wavefunction [see also Eq. (|40|) 1, the prime denotes the reduced Green function, and E is 
the energy of the raP state ("LAL" = loop-after-loop). The double-logarithmic term A-E^al 
originating from (p5|) reads 



AE- 



21og 
LAL 



a 
9 



In 2 



(Za) 2 m 



X ( <Pn, 



l,m 



P(H-E)P 

m m 



H-E 



P(H-E)P 

m m 



(26) 



In order to obtain this result, the denominator of the Green function H—E+iv has been expanded 
in powers of H — E within the integration region uj € [(Za) 2 m, e\. Using the commutator relation 



ABA = - ([A, [B, A]} + A 2 B + BA 2 ) 



(27) 



7 



with A = p/m and B = H — E, the matrix element can be rewritten in a much simpler fashion, 
and the double-logarithmic term becomes 



AR 



21og 
LAL 



I (-Y In 



9 Vvr 



{Za) 2 m 



1 

rn 



We have 



l> i, m \p 2 (H -E) p 2 \(f>nx 



, \X,l,m\p 2 {H - E)p 2 \(j} njlt 



(Zafm 5 (A 8 



5 15n 2 



(28) 



(29) 



Note that for S states, the above matrix element is divergent, and a regularization of the matrix 
element gives rise to an additional (triple) logarithm B^. With the natural ultraviolet cutoff 
e ~ m for nonrelativistic QED, we obtain from ( p8[ ) and (29) the following double-logarithmic 
contribution, 



A£g5(n,i = l) 



a\ 2 (Za) e m 

IT J 



TV' 



-] ^-ln 2 [(Z«)- 2 ] (1 



135n 2 



(30) 



Note that the presence of an additional double-logarithmic term originating from the loop-after- 
loop diagram in Fig. |l] (c) in the Coulomb gauge does not imply that the result given in (24) for 
the total value of Bq2 is necessarily incomplete, but it means that additional double logarithms 
have to expected if, e.g., this diagram is treated numerically, and numerical and analytic results 
are compared. For S states, an additional contribution to the triple logarithm B^ originating 
from the loop-after- loop diagram was found in [15,18, 19], but the result originally obtained in || 
for the total value of -E>63 was confirmed in [17, 22 1. In the following section, we will derive the 
result ( ]24| ) by an independent calculation which includes the entire gauge-invariant set of the 
diagrams in Fig. Q in a rigorous way. 



8 



5 Derivation Based on NRQED 



We start from the expression [see Eq. (16) of [22]], 

A-Enrqed = - ( » 2a o ) / dwi wi / duj 2 ^2 
\3irm A J Jo J 



F H — E + u>i H - E + uJi+uj2 H - E + uj 2 
1 / i 1 i 1 1 

+ - (p — p> — — ; p> 



2 \ H — E + uj\ H - E + uji+uj 2 H - E + uj 2 

1 / i 1 , 1 1 
+- — — p* — — p J 



2 \ H - E + uj 2 H - E + UJX+UJ2 H - E + u\ 




i 1 1 j 

" "' x P H-E + ui H -E + uj 2 P 

-^h(p' ' A — ^(p< ' A)- pi) 

LOl + U) 2 \ H — E + U 2 j U1+LU2 \ H - E + L0± /J 

All of the matrix elements are evaluated on the reference state \(f>), which can be taken as the 
Schrodinger wave function. 

Within the e-method ||[l2|, [24j], we extract those divergent contributions from (31) that in- 
volve double logarithms a 2 (Za) 6 ln 2 [e/(Za) 2 m\ (we may put e = t\ = t2 for simplicity). 
These logarithms correspond to the ultraviolet divergence of NRQED and are generated by 
the contributions of two infrared photons < e,uJ2 < e). The divergences in e necessarily 
cancel at the end of the calculation due to contributions proportional to ln(m/e) ln[e/(Za) 2 m] 
which are generated by intermediate integration regions {u\ > e,W2 < e), and by terms propor- 
tional to ln 2 (m/e) originating from high-energy virtual photons {u\ > e,0J2 > e )- The latter 
terms correspond to the infrared divergent terms proportional to ln 2 (A/m) of the electron 
form factors. For a discussion of the related cancellations in the context of the e method, 



we refer to |12|] and the Appendix of For the double logarithms, the dependence on e 
cancels between the low-energy, the intermediate and the high-energy regions according to 
ln 2 [e/{Za) 2 m] +2m(m/e) ln[e/{Za) 2 m] +ln 2 (m/e) =ln[(ZoT 2 ]. 



There are nine terms in curly brackets on the right-hand side of Eq. (31) which we would like to 
denote by T\ —Tg. These fall quite naturally into six groups, giving rise to six double logarithms 
L\ - Cq according to the following correspondence: 

• Tx -> Ci, 

• T 2 + T 3 -> £ 2 , 

• T 4 -> £ 3 , 
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• T 7 -> £5, 

• T 8 + T 9 -> £ 6 . 



After an integration in the logarithmic region oj\ G [(Za) 2 m, e\] and W2 € [(Za) 2 m, £2], the 
logarithmic contributions can be expressed by matrix elements, evaluated on the reference state, 
according to the following formulas (again, we put for simplicity e = e\ = 62): 



£1 



£2 



*] In- 

7T> 



^rin 2 

IT. 



In 2 




4(p 4 (H-E)p l p 2 ) 



2 (p» p 7 (if - E) p> p l ) - A{p l (H - E) p i p 2 ^ 
9m 1 



(Zaf_ 

£4 oc (p i (H-E)p i } = 0, 



(p 2 (H-E)p 2 ) 
9 m 4 



£5 



£ 6 



In 2 



In 2 




4(p* (H-Efp^ 



-A(p l (H — E) 2 p ix 
9 m? 



(32) 

(33) 

(34) 
(35) 
(36) 

(37) 



All of these matrix elements are finite when evaluated on P states and on states with higher 
angular momenta. In deriving these results, use is made of the integrals 1\ - I3 listed in Ap- 
pendix [B[ In particular, 1% is used in deriving £1, I2 is used in deriving £2, and Cq can be 
derived using I3. The double logarithm £3 corresponds to Eq. (p8|). Summing all contributions 
£1 - Cq, we obtain 



£ 



a 



In 2 



-) In' 

7T / 




2 (p* p 7 ' (if - E) p> p l ) - {p 2 (H -E)p 2 ) 



9 m 4 



in agreement with formulas (|T^) and (|l~8|). Here, use is made of the equality 

(p i p j (H - E) p> p*) = vr (Za) (A5^ (r)> + - (p 2 (H - E) p 2 ) , 



(38) 



(39) 



which is valid for P states and states with higher angular momenta and can be derived using 
the commutator relation (p7|). We thereby confirm that the additional double logarithm £3 
generated by the loop-after-loop diagram Fig. |l| (c) is cancelled by an additional contribution 
from £2 according to Eqs. d33| ) and (|39"1). 

As a byproduct of the current investigation, we obtain the rigorous result that B§2 vanishes 
for states with higher angular momenta / > 2. This is because the expectation value of the 
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effective potential (fl7|), when evaluated on hydrogenic D, F, G,. . . states, vanishes: states with 
higher angular momenta behave as r l for small r, where I is the angular momentum. We thereby 
confirm a statement made in || [following Eq. (5) ibid.] where it was pointed out that a formula 
analogous to (18) holds for all states with I ^ 1 [see the text following Eq. (5) ibid.]. 



6 Results and Conclusions 



The results of the current investigation can be summarized as follows: In Sec. |2|, we attempt 
to clarify the derivation and physical origin of effective potentials || used for the approximate 
description of self-energy corrections in leading logarithmic accuracy, and to provide a more de- 
tailed derivation of known double-logarithmic corrections to the Lamb shift. In Sec. ||, restrict- 
ing the discussion to P states and states with higher angular momenta, we rederive, within the 
effective-potential approach, known results || for the leading spin-independent double logarithm 
for P states as given by the Bq2 coefficient [see Eq. (p4|)]. In Sec. |], we show that nonvanishing 
double logarithms have to be expected from the loop-after-loop diagram if this non-gauge in- 
variant term is treated separately (e.g., within a numerical evaluation). By contrast, within the 
effective-potential approach, the double logarithm for this diagram vanishes (see the entry in 
column 2, row 4 of Table 1 of Q). In Sec. ||, we show that a rigorous derivation of Bq2 based on 
the entire gauge-invariant set of diagrams in Fig. [l] confirms the result (|24]) for the total value of 
Bq2- In particular, the additional double logarithm originating from the loop-after- loop diagram 
cancels when the contributions of all diagrams are added, and B^2 vanishes for all states with 
angular momenta I > 1. 

A reliable understanding of the problematic two-loop corrections is important for the determi- 
nation of fundamental constants from precision spectroscopy p5|] . We would also like to stress 
that analytic calculations, even in the low-Z region, could be supplemented by accurate numer- 
ical evaluations in the near future. Recently, a complete evaluation of the two-loop self-energy 
effect for high-Z has been reported [26|| . A comparison of the numerical to the analytic results 
represents a crucial test for both methods p7| . In order to provide for a reliable comparison of 
numerical vs. analytic results, it is helpful to thoroughly analyze and understand the logarithmic 
terms from each one of the diagrams in Fig. |l]. As outlined in Sec. 5 of ||, the most accurate 
theoretical predictions for the energy levels can be obtained using a combination of analytic and 
numerical results. 
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A Appendix A: Analytic Evaluation of a Matrix Element 



In this appendix, we discuss the derivation of the expression (19) 

|2~ 



A 



■ > n,l=l,m 



r=0 



for hydrogenic P states. In Eq. (18), the Laplacian operator acts on nonrelativistic, Schrodinger 
wavefunctions, which are given by 



4>n,l=l,m( r ) = Rnli^Yx 



(40) 
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where R n \(r) is the radial component, Y\ m {0,4>) is the spherical harmonics with the polar- 
coordinates r, 9 and <f> and with quantum numbers (n, I = l,m). Since the quantum number 
/ = 1 than the magnetic quantum number can be m = and m = 1,-1. For the sake of 
simplicity we consider the m = case, 



Y-i 



l m =0l<7, ' 



4vr J 



1/2 



cos ( 



The Laplacian in (18) can be written in polar-coordinates as 

1 d 



A A = | — — r 2 — ] + — 



■a 9 \ 1 32 



(41) 



(42) 



where A r corresponds to the radial component and Ag j, stands for the angular-dependent part 
of the Laplacian operator. One easily obtains 



A 



VZ=l,m=0 



(r)| 



AM 



4tt 



cos 2 6 + R 2 nl 



3_ 2_ 

47T r 2 



,1 — 3 cos 



(43) 



The final result ( |i9D should be independent of the angle 6, i.e. independent of the spatial direction 
in which the origin is approached, and independent of the magnetic quantum number. Therefore, 
we may postulate that the ^-dependent terms in (|4^) mutually cancel. Alternatively, we observe 
that since (|l9|) should be independent of the angle 6, so that so that we may set = it/ 2. 
Reading off the ^-independent part of (43), the following result can be obtained: 



A 



W,oO)l 



+ o 4tt lr 2 nl 



(44) 



The radial component of the Schrodinger wavefunction for hydrogenlike P states (R n i) is defined 
by the associated Laguerre polynomials (L^ +1 ) which read 



Rni(r) 



(n-2)! 



(n+ l)! 3 (2n)! 



1/2 



2 \ 5/2 

J r exp 



na B ) 



n+l 



2r \ 



L 3 n+1 (p) 



n+l 



j=0 



(45) 



where the Bohr radius is as = l/(Zam). Using this relation, it is straightforward to obtain 



A 



I0n,i,o(r)r 



+o 3tt 



[(Za) 



5 m 5 " 



n 



1 



n" 



(46) 



which is equivalent to Eq. (|T^). 



B Appendix B: Double— Logarithmic Integrals 

In this Appendix, we provide the results for certain integrals which may be used in order to 
extract the double-logarithmic contributions of order [Zoif ln 2 [e/(Za) 2 ] from the NRQED two- 



loop self-energy (31). We have two photon energies uj\ and u>2 an d denote arbitrary matrix 
elements of the various occurrences of the operator H — E, scaled by (Za) 2 , by the symbols 



12 



Ai, A2 and A3, respectively. The symbol ~ in this Appendix is meant to indicate that only the 
double-logarithmic terms of order (Za) e are selected. We have, 

f 62 , (Za) 2 1 {Za) 2 

dwi uj\ I auo 2 0J2 — — ; — A /r7 _\<) — — — ; — A lr7 ^\o ~ — ; — a iry^\ 2 



J 



u)\ + A x (Za) 2 w\ + uj 2 + A 2 (Za) 2 co 2 + A 3 (Za) 



i( Za)6u i^ ln (£? (Al+A3) ' (47> 

f (Zaf 1 (Zaf 



Jo + A\ (Za) 2 lji + uj 2 + A 2 (Za) 2 u> 2 + A 3 (Za) 2 

~ ^ Za)6ln w kl (^?^ + ' 43 - A2) ' (48> 

« /•« 1 (Za) 2 

dcji cji / auj 2 oj 2 



JO 



uji+uj 2 uj 2 + A (Za) 2 
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